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Non-adiabatic charge pumping through a single-level quantum dot with periodically modulated 
parameters is studied theoretically. By means of a quantum-master-equation approach the full 
counting statistics of the system is obtained. We find a trinomial-probability distribution of the 
charge transfer, which adequately describes the reversal of the pumping current by sweeping the 
driving frequency. Further, we derive equations of motion for current and noise, and solve those 
numerically for two different driving schemes. Both show interesting features which can be fully 
analyzed due to the simple and generic model studied. 
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I. INTRODUCTION 

Pumping of electrons through nano devices by a time- 
dependent modulation of device parameters has received 
a lot of interest over the past years. Such pumps are in¬ 
teresting for many app lications, but they are particularly 
useful in metrology 1 In this context, experimental and 
technologicalprogress has lead to high-frequency (in the 
GHz regime (P® an d high-accuracy charge pumping. 

Moreover, pumping is also interesting for addressing 
fundamental questions connected with the transport of 
quantum particles. Knowledge about the full counting 
statistics (FCS) of the pumped electrons allows a detailed 
understanding of relaxation and quantum effectsP. FCS 
for Coulomb-blockade systems has been studied theoret¬ 
ically in the context of stationarjHliiU driverP-hLU an d 
nano-electromechanicap3®23 systems. In particular, un¬ 
derstanding the noise of the pumping current is relevant 
for high-accuracy pumping. Accordingly, the noise in 
different setups of driven devices has been investigated 
theoretically 1 ' 19 and measured, for example, in a charge 
pumjPSl. 

While adiabatic pumping is very well studiecP^^H, the 
description of non-adiabatic effects remains challenging. 
In view of the experimental developments towards higher 
frequencies, this regime becomes increasingly relevant 27 . 
More over, non-adiabatic driving can lead to interesting 
effect^ 5 like the possibility to reverse the pumping 
current by sweeping the driving frequency^®. In order to 
exploit such effects a better understanding of the FCS 
for fast pumping is necessary. 

In this article we consider an electron pump modeled 
by a single-level quantum dot at zero bias. Based on the 
formalism given befor^f^ we calculate the cumulant gen¬ 
erating function at zero temperature and for arbitrary 
driving schemes. We show that in the adiabatic regime, 
the probability distribution for pumping electrons is al¬ 




ways binomial, which implies that the pumping is unidi¬ 
rectional. In the non-adiabatic regime we obtain a tri¬ 
nomial distribution, which also explains the occurrence 
of the current reversal®! We show that fluctuations at 
the reversal may show either a minimum or a maximum, 
depending on the driving scheme. 

The outline of the article is as follows. In the next sec¬ 
tion we derive the FCS for a quantum-dot electron pump. 
We also discuss the calculation of the pumping-current 
noise. In Sec. m we consider a specific time-dependence 
of the tunneling rates and the dot energy and compare 
the numerical results for this model to the analytic ex¬ 
pressions obtained in Sec. [TTJ Finally, we conclude with a 
summary and a discussion. 

II. THEORY 

We consider a resonant-level model® characterized by 
a time-dependent energy level e[t) and time-dependent 
couplings to the left and right reservoirs given by tun¬ 
neling rates Tl^) and r R (f), respectively. As indicated 
these parameters explicitly depend on time. The actual 
time dependence will be specified below. 

In order to obtain the FCS we consider the number 
of electrons N tunneling through the left barrier. Note, 
that in general this number is different for the right bar¬ 
rier. However, time-averaged quantities in the steady- 
state regime do not depend on the barrier (left or right) 
they are calculated for. The statistics of N can be found 
from the characteristic function <1>(x,f) = (exp[iy./V]), 
where x is the counting variable. Knowing <J>(%, t), one 
can deduce the moments ( N m {t )) by differentiation. 

Under the usual conditions®!, short reservoir 
correlation-time and Coulomb blockade, the char¬ 
acteristic function for the present model can be 
determined from the following equatiorP 


mm -/(t) [e-*r L (f)+r R (t)] \ , 

-/(*) [e ix r L (t)+r R (t)] mm ) ’ 
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where we have used 

^ ^ 1 + exp(e(f)/fc B T)’ 

r(t) = r L (t) + r R (t). 


f(t) = i-m, (ib) 

(lc) 


Note that the counting field x “measures” tunneling from 
and to the left reservoir only. For \ = 0; Eq. 0 reduces 
to the usual (Markovian) master equation for the diag¬ 
onal elements of the reduced density matrix 31 ! In this 
case, the components of the vector p = (pcnPi)* corre¬ 
spond to po(t) = 1 — n(t) and P\(t) = n(t) with n(t) the 
average occupation of the level at time t. The character¬ 
istic function is given in terms of the solution p of Eq. ([T]) 


$(x>0) = qp(x,</>) with q= (1,1). (2) 

The product is to be understood as a scalar product of 
two two-component vectors. 

Provided, that all external parameters change periodi¬ 
cally with frequency fl, we can substitute cj> = fit to gelf® 

d 1 

0 ^ p(x, </>) = -^M(/>)p(x, <j>) (3) 

with L x defined in Eq. (jTJ) . The solution after one cy¬ 
cle can be written as p(x, </> + 27r) = Ap (x,<P) with the 
matrix 


A = T exp 



where T denotes the time-order prescription. 

Assuming that the counting fields can be switched on 
and off adiabatically, the g ener ating function for k count¬ 
ing cycles is obtained froirPES) 

<f> fc = qA fc p, (5) 

where p is the steady-state solution for x = 0- For a 
large number of cycles k the char acter istic function is 
determined by the largest eigenvalu e) 9 ! 14 ! of A. 


A. Current and Noise 


For many applications it is sufficient to know the first 
two moments, Q = (N) and AQ 2 = (AA 2 ), of the 
pumped charge. The moments are given in units of e 
and e 2 . In the following, we will first derive a set of 
equations, which allow the calculation of both moments 
for arbitrary frequencies and temperatures. Since we as¬ 
sumed the counting to occur in the left reservoir, all time- 
dependent quantities refer to that reservoir and the label 
is suppressed. 

As will be shown, it is more convenient to consider the 
time-derivatives, I(t) = d t (N) and S(t) = d t (AA 2 ). 
Starting from the definition of the current, one gets the 
following expression 


d d d 

/(«)- S W-q 5gra . 


<9L y 

= 

x—o d( lx) 


x=0 


= r L (t) [f(t) - n(t)\ , 


( 6 ) 


where we have used Eq. ([!]) and qL x= o = (0, 0). Similarly one finds 


s(t) = -\( N2 )~( N r 


= q 


d d 2 
dt d(ix)‘ 


- 2 1 


<9(ix) 


p 

= q 


X—0 L 


= -2T h (t)r(t) + r L (f) [f(t)n(t) + , 


d 


- 21 - 


d 


5 2 L, 


d{ ix) 2 ^(ix) 9{ix) d{ i%) 


- 2 x 


x=o 


( 7 ) 


where 


equations are 


r(t) = p[{t) - n(t) \p' 0 (t) +p[{t)] , 


p'At) = 


dPj{x,t) 

d(i X ) 


x=o 


for j = 0,1. 


(8a) 

(8b) 


-^n(f) = F(f) lf{t)-n{t)] , 

(9a) 

J^r(t) = ~ T (t)r{t) 

(9b) 


+ r L (t) - lf(t)-n{t)} 2 ] . 


To calculate / and S , the time evolution of the occupation 
n and the newly defined quantity r are needed. Both of 
them can be obtained from Eq. ([l]). The corresponding 


Equations 0 can be solved numerically in a straight¬ 
forward way. The moments Q and AQ 2 are found by 
integrating / and S over one period, respectively. For 
low frequencies and finite temperatures, one can obtain 
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analytical expressions for both quantities. This case is 
discussed in Sec. Ill Cl below. 


B. FCS at zero temperature 

In the following section we calculate the generating 
function at zero temperature, i.e., when the Fermi func¬ 


tion only attains the values 0 and 1. To this end we as¬ 
sume that for 0 < cf> < n the level is charged (/= 1 , /= 0 ) 
and for n < <j> < 2n it is de-charged (/=0, /= 1). In this 
case, the solutions to Eq. 0 at the end of the respective 
half-periods can be obtained analytically. Overall, we get 
for one period 


n(v 9 7r \ _ a f v w v nl with A (e ix 7L+7R)(e ix /?L+/?R) e ix 7 L +7R \ 

P(x> 2 tt) - A o (x)p( X ,0) with A 0 (x) = (1 - 7L - 7 R)(e ix ^ L +^) 1~7 l-7r ) 


and the abbreviations 


a ( 1 

^ = / n ^-n“ exp “n 


^V r(^)j , 7q = fjdt ^ exp ^ V rtf )j . 


( 10 a) 


( 10 b) 


Note that P a and are completely determined by the 
time-dependence of the rates T a . They define proba¬ 
bilities for charging or cle-charging the level either from 
the left or the right reservoir, respectively. Obviously 
Pa, 7« > 0. Further we note that 


1 for x = 0 only to 


expression in Eq. (14) yields 
first order in p, 7 . 

Equation (14) characterizes a trinomial probability dis¬ 
tribution. Accordingly, there are three relevant probabil¬ 
ities 


Pl + Pr = 1 - exp J <1. ( 11 ) 

Since both quantities, /3 l and /3 r, are positive it follows 
from (11) that /3 l, Pr < 1. Similar arguments hold for 7 l 


and tr. In order to simplify the notation, in the following 
we use the definitions 


p = 1 - p L — /3 r , 7 = 1 - tl - 7R • (12) 


Finally, the eigenvalues of the zero-temperature matrix 
Ao(x)> which is dehned in Eq. © 0 , are given by 


A± = - 


six) ± \Zg 2 (x)-^P r r 


g(x) = P+ 7 + (e 1X 7 l +7r ) (e IX /?l +/3r ) 


(13a) 

(13b) 


from which the current and noise can be obtained by 
differentiation. The eigenvalue with the largest absolute 
value is A + . This expression is exact under the assump¬ 
tions stated above and can be used to obtain the FCS of 
the pumped charge. 

To gain further insight about the nature of the statis¬ 
tics, it is useful to consider limiting cases. If /?, 7 -C P, 7 
one can write down the following generating function 

($ fe ) 1 /fe ~7L/3Re- ix +/?L7Re ix 

+ (P + 7 + /?l7l + /3r7r) • (14) 


This limit describes the situation where the level is al¬ 
most fully charged in the first half cycle and correspond¬ 
ingly de-charged in the second half cycle. Note that the 


P- = 7l/3r, P+ = /3l7r and p 0 = l - P- -p+, (15) 

which describe a process of an electron being transferred 
to the left reservoir, to the right reservoir or no transfer, 
respectively. It follows that the average charge and noise 
per period are 

Q=P+~P~, (16a) 

AQ 2 = p + (l-p + )+p_{l-p_)+2p + p_ . (16b) 

For p + = p_ one finds that the pumped charge vanishes 
Q = 0, while the noise remains finite A Q 2 = 2 p + for 
p_l_>0. At th e vicinity of this point one obtains a current 
reversalAlAil 

In the opposite limit, P, 7 « 1, where the energy level 
is nearly empty during one cycle, one can approximate 
the generating function by 

($ fe ) 1/fc = c+ \/(e- ix 7L+7R)(e ix /3L+/3R) ■ (17) 

Here, c is independent of x and guarantees <f>fc| x= o = 1. 
Equation ( |17[ ) describes in general a complicated prob¬ 
ability distribution. However, if one barrier is dominat¬ 
ing the charging and the other one dominates the de¬ 
charging, one obtains a binomial distribution for a half¬ 
charge transfer. For example, if /3 R ss y L « 0 one gets 

(470 1/fc «c' + e ix/2 /3 L 7R. (18) 

This “fractional” behavior has been discussed previ¬ 
ously in the context of a nano-electromechanical charge 
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shuttle^. It reflects the fact that the cycles are no longer 
independent of each other. For many counting cycles, 
however, the behavior can effectively be described by the 
independent transfer of fractional charges. For later ref¬ 
erence we note that the charge per cycle can be obtained 
from Eq. © to read 


/3l7r ~ /3r7l 
2\/(/3l+/3r)(7l+7r) 


(19) 


C. Noise for low frequencies and finite 
temperatures 


At low frequencies, it is sufficient to consider the in¬ 
stantaneous contribution from Eqs. as done similarly 
before 2 ®. Consequently, the time-derivatives are set to 
zero and one obtains 


n(t) = f(t ), r (t)r(t) = T L (t)f(t)f(t) . ( 20 ) 


The instantaneous current and noise become 


m = o, 


S{t) = 


2 mm 


i/r L (t) + i/r R (t) ' 


( 21 ) 


Not surprisingly, the current is zero, since the level is at 
all times in equilibrium with the reservoirs. The noise 
is nonzero and is called equilibrium, or Nyquist-Johnson 
nois^H. It is caused by the statistical nature of the elec¬ 
tron occupation in the reservoirs characterized by the 
Fermi distribution. Using the property // = —k B T^ 
the expression for the noise becomes 


S(t) 


2 fcBT i/r L (f) + i/r R (f) de ' 


( 22 ) 


The noise per period is obtained from integrating over 
one period, 


A Q 2 = -2 


k B T 

n 



[de{(t>) / d(j>\ 1 [df(<j))/d<j)\ 

i/r L (</>) + i/r R (</>) 


(23) 


The integral gives a value, which is independent of f1 but 
depends on the details of the external driving, i. e., the 
time-dependent level and tunneling rates. Consequently, 
the noise per period will be proportional to k B T/kl. 


III. RESULTS 

We study the first two cycle-averaged moments of the 
FCS: the charge Q and its fluctuations A Q 2 for an elec¬ 
tron pump which shows current reversal and rectification 
effect in the non-adiabatic regim^H. They are obtained 
by integrating I[t) and 5(f), as given by Eqs. ([ 6 ]) and 
0 - respectively, over one cycle under steady-state con¬ 
ditions, i. e., p(<j>) = Ap {(j>) with A defined in Eq. Q. 


A. Exponential modulation 

To this end, Eq. 0 has been integrated numerically 
with the following pumping parameters 

e(t) = cos(flt), (24a) 

r L (f) = [r /2] exp ( 6 [cos (fit — <5) — 1]) , (24b) 

r R (t) = [r /2] exp ( 6 [cos (fit) — 1]). (24c) 

The same parameters have been used before 29 . They 
allow for the exponential dependence of tunneling r ates 
due to oscillatory gate-voltages as used in experiments 5 ®. 
The pumping is characterized by two parameters, the 
driving frequency fl and the time or phase delay <5 of 
the left barrier with respect to the right one. The latter 
one is locked to the level e which oscillates around the 
Fermi energy of both reservoirs /il = /i R = 0. If not 
specified differently, the temperature in the contacts is 

k B T = r/ioo. 

Figure [l] gives an overview of the results. We have var¬ 
ied the frequency fl in a rather large range (10 -5 ... 10 T) 
in order to cover the asymptotic adiabatic and non- 
adiabatic behavior for low and high frequencies, respec¬ 
tively. One can clearly see (Fig.[l£i) that the direction of 
the charge transport does not only depend on the phase 
delay S which controls the oscillation of the left barrier 



phase delay 8 


FIG. 1: Pumped charge Q and fluctuations AQ 2 as a function 
of phase delay S and driving frequency fl as obtained with the 
set of Eqs. 0 0 for the time dependence of quantum-dot 

level e and couplings r R>L given in Eqs. ([24]) . 
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FIG. 2: Charge fluctuations 0 obtained from the rate equations 0 for three selected phase delays 5 and fcBT=r/100. They 
are compared to thos e o btained for a trinomial probability distribution (161 with p_ = 7l/3r, p+ = /3l7r and a binomial 
distribution function (251 with p = \p+— p~\. 





FIG. 3: Charge fluctuations as shown in Fig.[2] here with a double-logarithmic plot and for two additional temperatures 
IcbT = r/1000 and ksT = T/IO, respectively. The dashed lines show thermal noise according to Eq. (231. 


with respect to the right one. Rather the frequency is 
equally important and may be used to cause a current 
reversal as was discussed in detail before^. 

The charge fluctuations (Fig.JljD) show a similar pat¬ 
tern as the average charge, albeit somehow inverted. Not 
surprisingly they vanish (white areas in Fig.[T o), when 
the charge transport is quantized either to the left (blue 
areas marked with A in Fig. 0 ) or to the right (red ar¬ 
eas marked with B in Fig. 0 ) direction. In regions, how¬ 
ever, where this is not the case the fluctuations are finite, 
reaching sometimes a value of A Q 2 = 1 / 2 . 

One maybe tempted to explain the charge fluctuations 
with the well-known expression 33 


A Q 2 = p(l-p), 


(25) 


with p the probability to transfer one charge per cycle. It 
is p = |<5| = \p+— P-\ with p± defined in Eq. (15). As we 
will see this binomial description is only applicable for 
finite values of Q. It fails in regions where the current 
changes direction, i. e., where Q= 0. There are two excep¬ 
tions where Q and A Q 2 vanishes simultaneously. One is 
the region marked with C in Fig.JT] and the other one is 


the high-frequency regime where Q 3 > T. We will discuss 
all regions in the following. To be more quantitative we 
have plotted the charge fluctuations for selected phase 
delays of <5=— 37t/4, 5=0 and 5=+7t/2, respectively, in 

Fig.0 

For S=— 37 r /4 there is a current reversal at fl ~ 
5xl0 _3 r (cf. inset in Fig.0). Along with the result from 
the rate equations (Fig. [Tjand thick gray line in Fig. [2]) we 
show the binomial expression (25) and the trinomial ex¬ 
pression from Eq. ©>)■ The values of p + and p- in the 
latter case have been calculated numerically. Note that 
this is an approximation since we have assumed k^T=0 
in defining p±. Both shot-noise models agree qualita¬ 
tively with the numerical result. However, at the current 
reversal (Q=0) charge fluctuations A Q 2 are on the order 
of 1/10, i. e., they do not vanish as predicted by Eq. (251. 


The trinomial description (|16| o) captures this behavior 
even quantitatively as Fig.[2|r shows. Thus, the fluctua¬ 
tions reveal that the vanishing net transfer is not due to 
the fact that there is no transfer at all. Rather transfer 
to the right and to the left do cancel each other exactly. 

The discrepancy between binomial and trinomial de¬ 
scription becomes even more obvious for 5=0 which is 
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FIG. 4: Charge fluctuations 0 for harmonic modulation ( |26| ) of the tunneling rates obtained from the rate equations 0 for 
thre e selected phase delays 5 and pBT=r/1000. They are compared to tho se obtained for a trinomial probability distribution 
(161 with p- = 7l/3r, p+ = /3 l7r and a binomial distribution function (251 with p = \p+— p~\. 


shown in Fig. 0- Due to symmetry ( p + = p_ ) the charge 
Q vanishes for all frequencies 52 and so does the binomial 
expression (25]). In contrast, the trinomial expression, 
which simplifies to A Q 2 = 2p + = 2p_, describes the 
fluctuations for all frequencies even quantitatively. Com¬ 
pared to the frequency-driven current reversal in Fig. 0 
the fluctuations are even larger here. 

The case 5=+n/2, shown in Fig.0, is characterized 
by a quantized and thus fluctuation-free transfer over a 
wide range of frequencies (52 « 10~ 5 ... 10 _1 r). 

It remains to discuss the case where Q and AQ 2 van¬ 
ish simultaneously. For large frequencies this is rather 
obvious, p + and p_ become small for 52 T since the 
transfer to either side can be neglected. It is less clear for 
the region marked with C. It can be traced back to the 
coupling to right contact: both /3 r and 7 r, as defined in 
Eq. (10 o), do vanish. If the level can neither be charged 
from the right (/3 r) nor de-charged to the right ( 7 r) the 
probabilities of transfer from ( p + ) and to (p_) the right 
contact do vanish as well. 

Careful inspection shows that for all three cases the 
fluctuations increase above the values predicted by the 
two shot-noise models. The reason is that thermal or 
Nyquist-Johnson noise AQ 2 herm starts to become larger 
than the shot noise. This is shown in Fig.[3]where thermal 
contribution according to Eq. (23) is shown separately 
with dashed lines. Clearly this occurs for larger temper¬ 
atures k^T at higher frequencies 52. The frequency, at 
which this takeover occurs, depends on the phase delay 
<5, as the three panels of Fig. [3] show. 


B. Harmonic modulation 


Taking, instead of pi]]) and (24 


T a (t) = r 0 + Ti cos(52f + 7 t/2 — 6 a ) 


(26) 


yields a harmonic variation of the coupling rates, which is 
frequently used for modeling electron pumps. The same 


setup has been investigated in view of a current reversal 
before 2 ^. Using the same approach as in the last section, 
we calculate Q and A Q 2 for To = Fi = F/20, Sr = 0 and 
Sl = 8. The results for a temperature of /cbT=F/ 1000 
are shown in Fig. [4] 

In all three cases, the fluctuations obtained from the 
trinomial probability distribution ( |16[) agree very well 
with the results of the rate equations 0. Comparing to 
Fig- E one sees that the results are qualitatively similar 
for i5 = 0 and S = tt/2. In the latter case and for har¬ 
monic modulation the fluctuations do not become zero 
in the adiabatic regime since the pumped charge is not 
quantized (cf. inset of Fig. |4}:) in contrast to the expo¬ 
nential driving (cf. inset of Fig. 0). 

The biggest qualitative difference is observed for S = 
— 37t/ 4. For harmonic modulation the charge fluctua¬ 
tions display a single maximum, which is attained in the 
vicinity of the current reversal. This is in contrast to 
the exponential modulation, where the fluctuations were 
found to have a minimum. This difference is connected 
with the behavior of the left/right probabilities p± in the 
vicinity of the reversal frequency 52*. It turns out that 
the slopes have the same magnitude but different signs. 
The curvatures may be different. Therefore we consider 
the following simple 52-dependence 

P± (fi) = P* ± p'* [n-n*] + p'*'± [52 — 52*] 2 (27) 


with p* = p±(52») an d Q (52*) = 0. For such a para- 
metrization the noise (16b) becomes up to 2nd order 


AQ 2 (52) « 2p* - [4p* 2 —p*][52 — 52*] 2 , 


(28) 


where the abbreviation p" = p" + +p"_ was used. The 
noise has an extremum at 52*, which can be a minimum 
or maximum depending on the relation of 4p* 2 and p*. 
Whereas in the harmonic case the curvatures of p± at 
52* have different signs but similar magnitudes, and thus 
P* ~ 0> the situation in the exponential case is differ¬ 
ent. Here one finds a non-linear increase/decrease of p± 
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with p” 4 p'J 1 . Correspondingly, close to the reversal 
we observe a maximum and a minimum of AQ 2 (fl), re¬ 
spectively. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have investigated the counting statis¬ 
tics for non-adiabatic pumping of electrons through a 
single-level quantum dot. For zero temperature we de¬ 
rived an analytic expression for the generating function 
[Eq. @] in terms of the probabilities for charging or 
de-charging the level during one pump cycle. In the case 
where those probabilities are large (the level is almost 
completely filled and emptied), we found a trinomial 
probability distribution for the charge transfer. The as¬ 
sociated elementary processes correspond to an electron 
being transferred to the left reservoir, to the right reser¬ 


voir or no transfer. This has the important consequence 
that the transferred charge per cycle can vanish while the 
charge fluctuations remain finite. It also shows that the 
current-reversal does not rely on interference effects. 

Our findings are corroborated by numerical simula¬ 
tions of the first two moments Q and AQ 2 for two driving 
schemes (exponential and harmonic). To this end we de¬ 
rived a set of ordinary differential equations, valid for 
arbitrary time-dependencies, which were solved numeri¬ 
cally. Those equations may also be used in connection 
with pulse-shaping techniques, which allow for optimiz¬ 
ing the pumping accuracjEl. 

Our calculations show that the driving frequency and 
the phase delay are important parameters, which both in¬ 
fluence the statistics of the pumping charge. This demon¬ 
strates that the ability to control the phase delays poten¬ 
tially provides an additional knob to improve the perfor¬ 
mance of electron pumps. 
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